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Abstract: We define and study the properties of generalized beam functions (BFs) and 
fragmenting jet functions (FJFs), which are fully- unintegrated parton distribution functions 
(PDFs) and fragmentation functions (FFs) for perturbative k±. We calculate at one loop 
the coefficients for matching them onto standard PDFs and FFs, correcting previous results 
for the BFs in the literature. Technical subtleties when measuring transverse momentum 
in dimensional regularization are clarified, and this enables us to renormalize in momentum 
space. Generalized BFs describe the distribution in the full four-momentum of a colliding 
parton taken out of an initial-state hadron, and therefore characterize the collinear initial-state 
radiation. We illustrate their importance through a factorization theorem for pp —> £ + l~ + 
jets, where the transverse momentum of the lepton pair is measured. Generalized FJFs are 
relevant for the analysis of semi-inclusive processes where the full momentum of a hadron, 
fragmenting from a jet with constrained invariant mass, is measured. Their significance 
is shown for the example of e + e~ — > dijet+/i, where the perpendicular momentum of the 
fragmenting hadron with respect to the thrust axis is measured. 
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1 Introduction 

In the description of high-energy scattering processes, the rigorous identification of the contri- 
butions from dynamics at different energy scales is achieved through factorization theorems. 
These allow one to systematically separate the short-distance behavior, which is calculable 
in perturbation theory, from process-independent, non-perturbative contributions. Further- 
more, factorization enables one to sum series of large logarithms of mass scale ratios, which 
would otherwise make the standard perturbative expansion unreliable. 

In inclusive processes with colliding hadrons, nonperturbative effects are encoded in par- 
ton distribution functions (PDFs). The standard PDF fi(x,fi) describes the distribution 
in momentum fraction x of a parton of type i = g,u,u,d, . . . inside an incoming energetic 
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r M = (— (x' — x)p ,—k±) 

Figure 1. Kinematics of the initial-state jet described by the generalized beam function. 

hadron. In this paper we will discuss fully-unintegrated PDFs, which depend on all com- 
ponents of the four-momentum k^ of the colliding parton. These allow one to keep exact 
kinematics rather than approximating incoming parton momenta as functions of x alone. 
Consequently, realistic distributions are obtained at leading order, for detailed final state 
measurements, like transverse momentum and invariant mass of jets [1]. The importance 
of fully-unintegrated parton densities was pointed out in the context of Monte-Carlo event 
generators in refs. [1-4], where they were called "doubly unintegrated parton densities" or 
"parton correlation functions". A field theoretic definition of fully-unintegrated PDFs and 
their applications in the pQCD formalism were discussed in refs. [5, 6]. 

Here we will focus on hadron-hadron processes where the colliding partons are far from 
threshold 1 and the initial-state radiation (ISR) emitted by a parton before entering the hard 
subprocess is constrained to a jet along the beam axis. This can for example be imposed 
through a veto on central jets or by an exclusive jet measurement, sec e.g. [7-9]. The momen- 
tum of the ISR is straightforwardly related to ft**, as illustrated in fig. 1 in terms of light-cone 
momentum components. 

The generalized beam functions B{(t,x,k±, fj,) are fully-unintegrated PDFs where the 
parton virtuality —t and transverse momentum k± are perturbative scales. They were defined 
in impact parameter space in refs. [10, 11] and named iBFs. They contain information on 
the initial-state jets concerning both their invariant mass and their momentum component 
perpendicular to the beam axis, which equals — k±. The real ISR pushes the transverse 
virtuality — t = k + k~ < of the colliding parton to be space-like 2 . In our kinematic 
setup {t, kj} <C Q 2 where Q is the hard scale associated with the partonic subprocess. The 
assumption {t, k 2 } 3> Aq CD allows us to calculate the dependence of the fully-unintegrated 
PDFs on t and k± in perturbation theory. The dependence on x can be written in terms of 

1 Here we do not resum logarithms of 1 — x, but these can be taken into account as shown in ref. [? ]. 

2 A proper definition of the beam functions requires a subtraction to remove the double counting of ultra-soft 
modes. From our private communication with the authors of refs. [10, 11] we learned that they were aware of 
this issue, but this was not addressed in their publications. 
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the standard PDFs, on which the fully-unintegrated PDFs are matched. In this paper we 
will consider t to be parametrically of the same size of k±, which avoids dealing with the 
resummation of logarithms in the ratio of these two scales. 

Generalized BFs extend the definition of the standard BFs Bi(t,x,fi) [7, 12], which 
have been extensively used to study processes with zero central jets [7, 13, 14], to include 
the measurement of the transverse recoil of ISR. Alternatively, starting from transverse- 
momentum dependent (TMD) PDFs (or unintegrated PDFs) [15? -22], one could think of 
the generalized BFs as a more differential version of them. However, the generalized BFs do 
not suffer from the rapidity divergences that affect TMD PDFs [23], as pointed out at 0{a s ) 
in ref. [10]. Instead, they have more in common with the standard BFs, like the anomalous 
dimension and one-loop matching onto PDFs, as we will show. 

One of our results is that by integrating a renormalized generalized BF over k±, one 
recovers the corresponding renormalized standard BF. No new UV divergences appear since 
the integration range of k± is bound by the measurement of t and x: 

t = —kk = xp = r^r > r, = k< , f 1.1) 

1 — x 1 — x 1 — x 

where k\ = —kj_k±^, > 0. Here is the total momentum of the ISR (see fig. 1), with r 2 > 
since this radiation is observed in the final state. Note that this relation between generalized 
and standard BFs is by definition true for bare quantities, but is nontrivial for renormalized 
quantities. By contrast, integrating the generalized BFs over t does not yield TMD parton 
densities. Here the range of t is not bound, so the integral over t generates new divergences. 

We study the renormalization of generalized BFs to all orders in perturbation theory. 
Their evolution is argued to be the same as for the standard beam and jet functions, which 
involves only the variable t. We discuss in detail how the proper definition of the relevant 
matrix elements is tied to the space where these get renormalized. In transverse momentum 
space, we show that in some cases a proper definition and renormalization requires us to mea- 
sure the transverse momentum k± in 2 rather than {d— 2) dimensions (all other momenta and 
polarizations are kept in d dimensions). We refer to the first case as CDR2 and the second 
case as CDR. In the transverse coordinate space both schemes are valid and lead to the same 
results for the generalized BFs. These impact-parameter-dependent BFs (iBFs) were intro- 
duced in ref. [10]. We stress that the evolution of the BFs for a perturbative renormalization 
scale does not change for nonperturbative t and k±, since the operator defining it is the same. 

In this paper we present the one-loop calculation of the matching coefficients between 
generalized BFs and standard PDFs, according to the hierarchy Q 3> k± ~ y/i S> Aqcd: 



Bi(t,x,kl,ii B ) = 2J / —r z v\ t ^—^ k l^B)fj{x ) ^B) 



-u,u,d,g 



A2 A2 

1 1 of QCD QCD ) 



(1.2) 



A matching of this kind was first discussed in the context of TMD fragmentation functions 
in ref. [24]. In refs. [7, 12] the analogous OPE was performed for the matching of standard 
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beam functions onto PDFs. In eq. (1.2) the matching scale [is should be chosen of the order 
of \ft to avoid large logarithms that spoil the convergence of standard perturbation theory. 
The Wilson coefficients Xij describe the (real and virtual) radiation building up the ISR jet 
from the parent parton j until the parton i enters the hard interaction, which has the same 
interpretation as for the matching of standard BFs onto PDFs [14, 25]. We work out Xij for 
i and j being either a quark or a gluon. We will argue that the remaining ingredients (hard 
and soft functions) in factorization theorems for zero central jets are the same as those used 
in conjunction with standard BFs. 

Our study of matching and renormalization allows us to compare with existing results 
obtained for iBFs in refs. [10, 11]. The diagonal coefficients X qq and X gg in eq. (1.2) agree with 
the ones obtained in ref. [11] but the mixing terms X qg and X gq do not. Taking the results 
of refs. [10, 11], and correcting the average over the incoming polarizations in dimensional 
regularization and fixing a sign resolve the discrepancy. 

In a fashion similar to beam functions, the fragmentation of a light hadron h within a jet 
originating from a parton i whose invariant mass s is constrained, is described by fragmenting 
jet functions (FJFs) Q^(s,z,fi) [26]. Here z denotes the hadron-parton momentum fraction. 
Fragmenting jet functions have the same infrared structure as the standard fragmentation 
functions (FFs) Dj(z,fi) and can be perturbatively matched onto the latter [27, 28]. In this 
framework, we have recently analyzed up to NNLL the cross section for e + e~ —> Xh on the 
T(4S) resonance where one restricts to the dijet limit by imposing a cut on the event shape 
variable thrust. This is relevant for the study of light-quark fragmentation performed by 
the Belle collaboration [29]. For large values of thrust we found that going beyond leading 
order in the cross section is important for a reliable extraction of the fragmentation function 
parameters [27]. 

Following our discussion on generalized BFs Bi(t, x, k±, fx), we also study the features of 
generalized FJFs Q^(s,z,p^-,fi), which also depend on the momentum pv~ of the observed 
hadron perpendicular to the jet axis. We discuss their renormalization and calculate the one- 
loop matching coefficients onto standard FFs according to the hierarchy Q 3> pr~ ~ yfs 3> 
-A-QCD- We argue that to all orders in perturbation theory the running of generalized FJFs is 
the same as that of the jet functions (as is the case for the standard FJFs). 

The paper is organized as follows. In sec. 2 we set up the theoretical framework, give 
the definitions of quark and gluon generalized BFs and FJFs, and discuss their all-order 
renormalization, both in momentum space and in impact parameter space. We make a 
general comparison with the existing SCET literature on generalized BFs [10] in sec. 2.2.3, 
and discuss the disagreement on the matching coefficients in sec. 3.2. In app. A we stress the 
need for ultra-soft zero-bin subtractions in the definition of BFs via a detailed calculation of 
the quark generalized BF where UV, IR and rapidity divergences are taken care of by different 
regulators. Sec. 3 is devoted to the calculation of the quark and gluon matching coefficients 
onto the standard PDFs/FFs. The details of the one-loop calculation of the gluon BF in 
momentum space using CDR2 are given in app. B. In sec. 4 (and sec. 2.2.4) we illustrate the 
relevance of generalized BFs and FJFs with factorization theorems for transverse momentum 
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dependent distributions. We discuss ^-distributions in Drell-Yan with a central jet veto and 
in gg — > H — > W + W~ — > l + vl'~v' + jets. We also analyse single-hadron fragmentation 
in e + e~ where a cut on thrust is used to restrict to dijet final states and the full hadron 
momentum is measured. 

2 Definition and Renormalization of Generalized BFs and FJFs 
2.1 SCET Ingredients 

As explained in the introduction, in this paper we deal with processes governed by three 
different perturbative scales: a hard scale Q associated with the partonic subprocess, an 
intermediate jet scale given by the square root of the jet invariant mass s or the transverse 
virtuality t, and a soft scale of order s/Q or t/Q, respectively. We will always consider the 
mass of the incoming/ fragmenting hadron to be negligible. We will also include effects at 
the non-perturbative scale Aqcd through standard PDFs and FFs. We restrict our analysis 
to the case that the transverse momentum of the ISR or the momentum of the fragmenting 
hadron perpendicular to the jet axis are of the order of the jet scale. After integrating out 
the hard dynamics, we are left with degrees of freedom describing the energetic radiation 
inside well-separated jets and the soft emission between them (collinear and ultra-soft modes, 
respectively). Therefore, Soft- Collinear Effective Theory (SCET) [30-33] is well suited to 
study generalized BFs and FJFs. 

Collinear modes are characterized by having large energy and small invariant mass, and 
therefore are conveniently described using light-cone coordinates. We introduce a light-cone 
vector whose spatial part is along the collimation axis, and another light-cone vector n M 
such that n 2 = n 2 = and n • n = 2. Any four- vector p^ can then be decomposed as 
= ,Pj) with p + = n ■ p, p~ = n ■ p and p^_, which contains the components of p^ 

perpendicular to and n^. The momentum of a collinear particle scales as {p + ,p~ ,Pj) ~ 
p~(A 2 ,l,A), where A <C 1 is the SCET expansion parameter. For the ultra-soft degrees of 
freedom, the momentum scales like = (q + , q~,q±) ~ p~(\ 2 , A 2 , A 2 ). 

The SCET fields for n-collinear quarks and gluons, £ n ,p e (y) an d A n ^ Pl (y) respectively, 
are labeled by the direction n and the large (discrete) momentum p£. Their argument y is 
conjugate to the small residual (continuous) momenta 3 . We explicitly exclude the case p^ = 
in collinear fields to avoid double-counting the soft degrees of freedom (which are described 
by separate ultra-soft quark and gluon fields). In practice, when calculating matrix elements, 
this is implemented via zero-bin subtractions [34]. 

Collinear operators are built out of products of fields and Wilson lines that are invariant 
under collinear gauge transformations [31, 32]. The basic building blocks are the collinear ly 
gauge-invariant quark and gluon fields, defined as 

Xn (y) = wl(y) Uv) , K±(y) = \ K(y) ^ n± w n ( y )] , (2.1) 

3 Throughout this paper, for notational convenience, we combine residual p r and label pi momenta into a 
continuous variable p: for example, 8^ - S(p~) = S(lj — p~) and similarly for the perpendicular components. 
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where iDf* , = , + gA*t , is the _L-collinear covariant derivative. The collinear Wilson line 



W n {y) 



perms 



V, 



n-A n (y) 



(2.2) 



sums up arbitrary emissions of n-collinear gluons which are 0(1) in the power counting. 

At leading order in the SCET power expansion, the interactions of ultra-soft gluons with 
collinear fields exponentiate to form eikonal Wilson lines. The ultra-soft gluons can thus be 
decoupled via the BPS field redefinition [33] 



X ^(y)=Y^y) Xn (y), 
B^\y) = Y^y)BUy)Y n (y). 



(2.3) 



The collinear fields we consider in this paper are those after this decoupling, and we drop 
the superscript (0) for notational convenience. Here Y n (y) is an ultra-soft Wilson line in the 
fundamental representation 



Y n (y) = Pexp 



,o 

/ du n-A us (y + un) 

J — oo 



(2.4) 



The symbol P in eq. (2.4) denotes path ordering of the color generators along the integration 
path. 

2.2 Generalized Beam Functions 
2.2.1 Definition 

The bare quark and gluon beam functions (BFs) are defined by the following proton matrix 
elements of bare SCET operators: 



B^%t,x,kj 



(2.5) 



v(i.--)(, > „(p-)\\„wi)(i -i- P +)i[s(k- 



1 



Vn)-5(kl-V* ± )Xn(0)] \ P n(p-)), 



7T 



a 



baic (t,x,k ± ) = 

- k- e(k-)(p n (p-)\B^ c ± (o) s(t - k-p + ) [s(k- - v n )5 2 (k± - v n ±) ^i(o)] \pn{p-)) 



where an average over the proton spin is assumed. The light-like vector is chosen such 
that the proton states have no perpendicular momentum, = p~n >1 /2. By boost invariance 
along the z-axis, these functions only depend on the momentum fraction x = k~ /p~, the 
transverse virtuality —t = k~k + of the parton, and the transverse momentum k_\_. The label 
momentum operators V act on the fields inside the square brackets whereas p + measures the 
plus momentum of any intermediate state. 

If we switch from momentum to position space, the fields in eq. (2.5) will be separated 
by a distance y M . Ultra-soft Wilson lines along the large y~ separation are factored out in 
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the BPS-redefined fields \n an d B n ± of eq. (2.5). The collinear Wilson lines in eq. (2.1) are 
along the small y + direction. Since y± 7^ 0, the BFs defined in eq. (2.5) are invariant under 
(ultra-soft and collinear) covariant gauge transformations because in this case the gauge field 
vanishes at infinity. The subtleties arising in the definition of these functions in singular gauges 
are similar to those affecting TMD PDFs, see e.g. refs. [23, 35? -38]. Since we deal with 
factorization theorems for gauge invariant quantities, we can always choose a non-singular 
gauge to carry out our calculations. 

The quark BF is a Lorentz scalar which only depends on the magnitude of k± and hence 
we could replace S 2 (k± — V n ±) with 5(k 2 — V 2 A/it. In the gluon BF, the measurement of 
k± allows a new Lorentz structure [10], namely 

B^(t, x,k ± ) = B^t, x, k 2 ± ) Lf + B 2 (t, x, k 2 ± ) L% u (k ± ) , 

Lf = \ , L^(k ± ) = ^ L + 9 f, L^L 2 ^k ± ) = . (2.6) 

Thus the tensor structure of the gluon BF and the dependence on the full k± must be kept. 
In a factorized cross section for two colliding gluons the Lorentz indices between two gluon 
BFs get contracted, see eq. (2.12). The new Lorentz structure L 2 only starts at one loop in 
the BF and it only starts to contribute to the cross section at two loops, since L\ ■ L 2 = 0. 

In order to interpret eq. (2.5) for the quark BF in QCD, in the Fourier transform, the 
field Xn(y) should be replaced by V(y,n)ip(y) with V(y, n) defined as Pexp ( — ig dsn • 
A(y + sn)) . In addition, soft subtractions need to be performed to avoid double counting of 
overlapping momentum regions [? ], in analogy to SCET zero-bin subtractions. 

An alternate definition of fully unintegrated parton distribution function, dubbed parton 
correlation function (PCF), was given in ref. [5] in the context of QCD factorization. In 
covariant gauges and in coordinate space, the PCF is defined as 

X, s _ WQ/)V%,n s )jy(0,n s )V>(0)|p) 

{V ' VP) Vs) (0\V(y, nr)Vt(y, n s )V(0, n a )Vt(0, n T )\Q) ^ <} 

which depends on all four components of the coordinate y, proton rapidity rj p and a soft "cut- 
off" rapidity r] s . n s = (—e v %e~' ns ,0±) is a space-like vector required to provide a boundary 
between left-moving and right-moving partons. ut = (— e - ' 2 *', 1, 0j_) is a non-light-like (or 
approximately light-like) vector and is required to regulate the rapidity divergences. Qualita- 
tively, the denominator serves the same purpose as zero-bin subtractions in SCET. However, 
there are important distinctions between our definition and the PCF. Note that the numer- 
ator in eqn. (2.7), which can be considered as a naive definition of the PCF, has space-like 
Wilson lines unlike our definition of BFs which has light-like Wilson lines. Another important 
distinction arises in the RG equations. Apart from the usual RG equation in u, the PCF 
satisfies a Collins-Soper equation in n s . This n s dependence would cancel against another 
collinear sector, so that the cross-section does not depend on rj s . In our case, the soft over- 
lap of the beam functions is only with momenta of order Q(A 2 , A 2 , A 2 ), hence BFs are not 



-7- 



rapidity divergent after zero-bin subtractions have been implemented. This is discussed in 
Appendix A. As a consequence, the standard RG equation in n is sufficient to sum all the 
large logarithms in our BFs. 

2.2.2 Treating Transverse Momenta in Dimensional Regularization 

In this section we discuss a technical issue arising in dimensional regularization when the full 
perpendicular momentum is measured rather than its norm. This is the case, for example, 
of the gluon generalized BF. We now show that the proper definition of this function in 
momentum space involves the measurement of k± in 2 dimensions (CDR2) rather than in 
d — 2 dimensions (CDR). With renormalization carried out in impact parameter space, CDR 
and CDR2 are both valid schemes which yield the same results for the calculations in this 
paper, as shown in sec. 3.2. However, we expect that CDR and CDR2 would in general lead 
to different results for the finite terms of other transverse-coordinate-dependent functions. 

The essential point of our argument is that the matrix element of operators (like the 
BFs) should be in integer dimensions, to have an unambiguous expansion in e = (4 — d)/2. 
This cleanly separates the divergences, facilitating the resummation of logarithms sector by 
sector. Requiring the matrix elements to be in integer dimensions fixes the dimension of the 
(5-function measuring the transverse momentum. 

We point out that the bare matrix elements have the same integer mass dimension as the 
renormalized ones. This follows from 

(O rcn ) = Z a l ® (O bare ) = (zjf Z% 2 • • • ) ^ ® (fafc ■■■), (2.8) 

where Z^. relates bare and renormalized ((pi) fields and Zq is the operator renormalization 
factor. The mass dimensions of the Z's and of the integration measure in the convolution 
"(8)" cancel each other, as is clear at tree level. 

Let us consider the example of eq. (2.5). For the fields, [% n ] = 3/2 — e and [B n ±] = 1 — e. 
For a single particle in d dimensions, [\p)\ = — 1 + e. In eq. (2.5) the e-dimensions thus 
cancel between fields and external one-particle states. Therefore the transverse momentum 
for the generalized gluon BF has to be measured in 2 dimensions, thus CDR2 is a sensible 
scheme. We stress that this would not be the case in CDR, where eq. (2.5) would involve 
a 5 d ~ 2 (k± — V n ±). For the quark generalized BF, which depends on the norm \k±\, both 
schemes will work. In CDR this requires using the identity 

S(ki - T^J/tt = r( ( 1 ^ ) ~j^ 5 2 ~ 2 *(k ± - V n ±) (2.9) 

In impact parameter space the gluon generalized BF is renormalizable both in CDR2 
and in CDR. In the latter case one performs a (d — 2)-dimensional Fourier transform with 
respect to k± [10, 11] which ensures that the bare iBF is in integer dimensions. The proper 
modifications of the Lorentz structures in eq. (2.6) for the bare beam functions in the CDR 
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and CDR2 schemes are 



CDR 





d- 



L^(k ± ) 




CDR 2 





d- 



L^(k ± ) 



(2.10) 



where k^_ is (d — 2) dimensional in CDR and two dimensional in CDR2. Here denotes 
the purely two-dimensional piece of g^. Since in CDR2 only the purely two-dimensional 
perpendicular momentum gets measured, the — 2e-dimensional contribution is not associated 
to a specific direction and should only appear in L\ . This is shown in the explicit calculation 
for the gluon BF in app. B. 

2.2.3 General Comparison with SCET Literature 

Generalized BFs are purely collinear matrix elements composed of collinear fields and states. 
Integrations over full phase space and loop momenta overlap with the momentum region 
corresponding to the ultra-soft modes. To attain a proper definition of collinear matrix 
elements we must implement zero-bin subtractions [34] . Previously, generalized BFs in impact 
parameter space (iBFs) were defined without explicit ultra-soft zero-bin subtractions [10, 11] 
4 . We show in app. A that in absence of these subtractions the generalized BFs suffer from 
rapidity divergences, by explicitly calculating the quark BF at one- loop with IR and rapidity 
regulators different from dimensional regularization. 

Generalized BFs naturally arise in the context of observables which are sensitive to 
collinear and ultra-soft modes. In our framework the initial-state radiation is constrained 
to energetic jets (described by collinear fields) with only ultra-soft radiation between jets 
due to our measurement of the transverse virtuality —t ~ Q 2 X 2 . Therefore the soft modes 
[which scale like Q(X, A, A)] cannot appear here as real radiation. Virtual exchange of soft glu- 
ons could in principle transfer perpendicular momentum between the two colliding partons. 
However, there is no measurement sensitive to this momentum transfer, rendering scaleless 
integrals for the soft contribution. In other words, a collinear particle within a jet can only 
recoil against another particle within the same jet. Our observables will become sensitive 
to soft radiation only if parametrically k 2 <C t ~ XQ 2 . However, in this case the contribu- 
tion of collinear radiation to t is power suppressed and one would have to deal with TMD 
PDFs rather than generalized beam functions. Alternatively, the generalized BFs do not have 
support on the soft region. We have assumed that effects from Glauber modes cancel out. 

By contrast in refs. [10, 11], the iBFs were used for cross sections where there is no 
constraint on t but only on k±. This in principle permits contributions from soft modes, 
which they treat as explicit degrees of freedom in the effective theory. In light of the previous 
paragraph this raises concern about the proper accounting of modes and power counting 

4 In our private communications the authors of these references agree that the proper definition of the iBF 
requires ultra-soft zero-bin subtractions. 
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in their factorization theorem. For example, their soft function depends on all momentum 
components, so it naively has an overlap with the ultra-soft region. Since the ultra-soft 
zero-bin subtractions were performed for the iBFs, which also depend on all momentum 
components, it is not clear to us how the ultra-soft subtractions for the soft function should 
be performed in their framework. A detailed discussion of their factorization theorem is 
beyond the scope of this paper. 

In sec. 3.2 we will resolve a discrepancy in the literature between refs. [11, 14] regarding the 
matching coefficients X qg and X gq of BFs onto PDFs, in favor of ref. [14]. In pure dimensional 
regularization zero-bins vanish and therefore the matching coefficients obtained in ref. [10, 11] 
should coincide with our calculation in this paper. As we will explain, the discrepancy is due 
to an oversight in ref. [11]. 

2.2.4 Renormalization 

We will now argue that the renormalization of the generalized BF is the same as that of 
the standard BF, to all orders in perturbation theory. In turn this equals the jet function 
renormalization [25] for which the anomalous dimension is known to three-loop order [14, 39, 
40]. This fact is quite useful since resummed calculations require anomalous dimensions at 
higher order in a s than the fixed-order contribution. 

Assuming the validity of a factorization theorem, we give an argument for the all-orders 
renormalization of the generalized BFs. In essence, in the factorization theorem the additional 
transverse momentum dependence only appears in the BF, and thus the generalized and 
standard BF have the same anomalous dimension. Consider, for example, the case of Higgs 
production through gluon fusion (gg — > H), where a central jet veto is imposed through 
the beam thrust event shape 7^ m defined in the hadronic center-of-mass frame [14]. This 
example involves the two different Lorentz structures in the gluon BF. It will be convenient 
to separately measure the contributions to 7~ cra from each of the two hemispheres orthogonal 
to the beam axis defined as 



T cm = T a + T b , T a = J2 e~ m , % = £ 0{- ni )\pt\ e* , (2.11) 



where the sum on % runs over all particles in the final state except the Higgs. Here \p^~ | and rji 
are the transverse momentum and pseudorapidity of the particle i with respect to the beam 
axis. The factorization theorem for small beam thrust in ref. [7] can be generalized to the 



case where we also measure the transverse momenta of the two ISR jets fc, 



a 



da 




dTadTbdk^dk^dY 




H J Bg (t a , x a , k a , /i) Bg ^(tb, Xb, , /i) 
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x I ^Bi(t a ,X a ,ka ± ,[i)B 1 (t b ,X b ,kb ± ,fl) + 



B 2 (t a ,x a ,k 2 ± ,fj,)B 2 (t b ,x b ,k b± ,fi)}, 



k 2 k 2 2 

K a± K b± 



where <jq is the Born cross section. The rapidity Y and mass rein of the Higgs are related to 
the momentum fractions x a ^ b by 

x a = ^-e Y , Xb = T ^-e~ Y . (2.13) 

"cm "cm 

The second equality in eq. (2.12) is obtained by inserting the Lorentz structures of eq. (2.6). 
The hard function H gg describes the virtual corrections at the hard scale mj/, and S^ emi is 
the incoming hemisphere soft function. Since only the large momentum components enter 
the hard function H gg , the additional measurement of the transverse momenta k a ± , k b ± does 
not affect it. By requiring that parametrically {T a ,T b } ~ {k a ±,k b ±}, the contribution of the 
ultra-soft radiation to the transverse momenta is power suppressed, so Sf^- is the same soft 
function as in ref. [7]. Soft degrees of freedom do not contribute, as explained in sec. 2.2.3. 

The //-dependence of the cross section cancels up to the order one is working, implying 
that the anomalous dimensions of the hard, soft and beam functions cancel each other. Since 
H gg and S^ emi are unchanged by the additional measurement of the transverse momenta, the 
products of the standard beam functions B g (t a , x a , y)B g (t b , x b , fx) and the products of the 
generalized beam functions Bi(t a , x a , k 2 , fi)Bi(t b , x b , k 2 , fj,), B 2 (t a , x a , k 2 , fj,)B 2 (t b , x b , k b 2 , y) 
all have the same running in fj>. The relative contribution of the B\B\ and B 2 B 2 Lorentz 
structures can be varied through the angle between k a and k b , so they are independent 
and cannot mix under renormalization. Since the variables t a ,t b ,x a ,x b ,k 2 ,k 2 can be all 
independently varied through T a ,T b ,mH,Y, k 2 , k 2 , it follows that the generalized gluon BF 
has the same renormalization as the standard gluon BF. The same is true for the quark beam 
function. 

The renormalization of the standard beam function equals that of the jet function [25], 
so we conclude that the renormalization of the generalized beam functions is given by 

B^ TC (t,x,k ± )= fdt' ZUt-t^^Biit^x,^^), Z i B (t,ri = Z i J (t,n), 
Jo 

V- r B i {t,x,k ± ,fi)= I d^nf%(t-l/,ii)Bi(lf,x,k ± ,ii), 7M*.A»)=7j(t,A»), ( 2 - 14 ) 



d/i jo 

where Zj and 7j are the jet function renormalization factor and anomalous dimension. There 
is no mixing between different parton types, so i on the right-hand side is fixed (not summed). 
Explicit expressions for 7^ may be found in app. D.2 of ref. [8] and for 7^ in app. B.3 of 
ref. [14]. 

It may be surprising that the renormalization in eq. (2.14) depends on i = —k~k + , 
instead of the Lorentz invariant combination t + k 2 = —k 2 . In SCET, Lorentz invariance 
is broken by the choice of and Instead we have reparametrization invariance (RPI) 
[41, 42], which encodes the arbitrariness in choosing and n^. RPI can be divided into three 
types: RPI-I and RPI-II transformations correspond to rotations of n and n, while for RPI-III 
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Figure 2. Fragmentation of the hadron h inside the jet Xh. 



n n g a n n ^ — )• e~ a . These transformations preserve the SCET power counting and 
the defining relations n 2 = n 2 = 0, n ■ n = 2. Let us now study the implications of RPI 
on the renormalization of the quark beam function. To this end we strip the beam function 
in eq. (2.5) of the external proton states, since the states do not affect the renormalization. 
Starting from the most general structure for the renormalization of the resulting operator O q , 
we find 



Og°(k~ ,k + ,kl) = dk-dk + dklZq(k-,k + ,kl,k-,k + ,kl,n)O q (k-,k + ,ki,ii) 



dk~dk + dk 2 



5(l - y;)s(k± ~ kl)Zf{k-k + ,k-k + ,fi) 

O q (k-,k + ,k 2 ,Li). (2.15) 



+ Zl^,k 2 ,k 2 ,fi 

The RPI transformations change V, V^x and p + inside O q into one other. By a suitable 
change of variables for k~ ,k + , k 2 and k~ , k + , k 2 , the operators take their old form and the 
effect of the RPI transformation is moved entirely into the Z factors. Based on RPI two 
structures are allowed: in Zf the renormalization is purely in t, in Z| the renormalization 
can depend on both k 2 and z. Eq. (2.14) shows that the structure Z| is absent to all orders, 
as verified by our one-loop calculation. 

2.3 Generalized Fragmenting Jet Functions 

The fragmenting jet functions (FJFs) Gi(s, z) describe the fragmentation of a light parton i to 
a light hadron h within a jet originating from i, where in addition to the momentum fraction 
z, the invariant mass s of the jet is measured [26]. Here, we will consider generalized FJFs, 
which depend also on the transverse momentum of the hadron ~ ^/s with respect to the 
jet axis, see fig. 2. We could define the generalized gluon FJF with uncontracted indices, as 
we did for the gluon BF in eq. (2.5). However, the new Lorentz structure only contributes 
at higher orders, especially for e + e~ collisions, due to the need of contracting indices with 
another gluon matrix element. We therefore restrict ourselves to the case with contracted 



Lorentz indices, for which only the magnitude \p^ \ is relevant. The transverse momentum 
measurement can be described by inserting in the definition of the standard in refs. [26-28] 
the additional function S 2 (p± — Ph±) = $(p± ~Ph±)/ 7T - 

The quark and gluon generalized FJFs are thus defined as 



gJ baxe ( S , z,pl) = 4(2vr) 3 [ d d Ph 5(pl) 5{p~ h - zk~) U{p 2 - p 2 ± ) £ 

J x c 

x tr[| (0\[S(k- - V n ) S d - 2 (V^ )xn(0)] S(k + - p+)\Xh){Xh\x n (0)\0) 
g* haie (s, z,pl) = -4(2vr) 3 k-J d d Ph 5(p 2 h ) 5{p~ h - zk~) U(pl - p 2 ± ) £ 



(d-2)(JV?-l) 



x(0|[5(A ; --P n )^- 2 (^)<i a (0)]5(A ; +-p+)|X/ t )<X/ l |^(0)|0), 

(2.16) 

where s = k 2 = k~k + is the invariant mass of the collinear radiation building up the jet 
(including the observed hadron). The label momentum operators V and act on the fields 
inside the square brackets, and the residual momentum operator p + acts on the intermediate 
states. The state \Xh) = \Xh(ph)) contains a hadron with momentum ph, and a sum over the 
polarizations of h is assumed. At variance with ref. [26], we choose to combine the label and 
residual minus and perpendicular momentum components into continuous variables (see also 
ref. [27]). In the trace is taken over color and Dirac indices, and the factor 1/(2N C ), where 
N c = 3 is the number of colors, comes from averaging over the color and spin of the parent 
parton. In Qg an average over colors and the (d — 2) polarizations of the gluon is performed. 

In ref. [26], a simple replacement rule was obtained that allows one to obtain the fac- 
torization theorem, when the momentum fraction z of a hadron in a jet is measured, from 
the corresponding inclusive case. Assuming cancellation of Glauber contributions, this can 
be extended to the situation where \p^ \ is also measured, 

Ji( s ) — ► 2 (2tt) 3 g ^ S ' Z, ^ A -'^ dz ^ d Ph± ■ ( 2 - 17 ) 

This equation is consistent with the fact that FJFs occur in SCETj factorization theorems 
where the transverse momentum of collinear particles is much larger than that of ultra-soft 
momenta. Soft modes do not contribute for the same reason as in the beam function case. 
Therefore the observed hadron only recoils against the other collinear particles within that 
jet. The replacement rule in eq. (2.17) still holds in the case of two (or more) gluon FJFs, 
if the angles of ph± are not observed. In this case the second Lorentz structure does not 
contribute. 

As a consequence of eq. (2.17), the jet function Jj and the FJF have the same renor- 
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malization, and thus the same anomalous dimension, 

£# baie (s, z,p* ± ) = fds' Z g (s - s', M ) Q*{sf, z,p* ± ,Li) , Z g (s, /i) = Z*j(8, fj.) , 
Jo 

Pty 3 * = f^ 8 ' 1 ^ 8 ~ s',»)gHs',z,Phx^) , 7M s 'M) = T}(s,M) , (2-18) 

where the parton type labelled by the index i is not summed over. In analogy with the 
argument for the all-order renormalization structure in sec. 2.2.4, let us take as example the 
decay rate of the single-inclusive jet-like process B — > {Xtt) u (.v, where the leading-power 
factorization theorem has the structure [26] 

dT = Hg*®S. (2.19) 

The measurement of the hadron _L-momentum affects only one function, namely Q£, the hard 
coefficient H and the soft function S do not change. Since the anomalous dimensions of all 
the functions of the RHS of eq. (2.19) must cancel, the _L-momentum of the hadron cannot 
change the anomalous dimension of Qf. More intuitively, the integral over p?, is bound once 
s is measured since p^ ± = p^p^ < k~k + = s. Therefore integrating over p£ ± cannot produce 
UV divergences. 



3 One-loop Matching 

In this section we discuss the matching of beam functions onto PDFs and fragmenting jet 
functions onto fragmentation functions. The one-loop matching coefficients are presented. 
We also provide the conversion relation from perpendicular-momentum to impact-parameter 
space. We compare to previous results for the matching coefficients in refs. [10, 11] and discuss 
the disagreement. 

3.1 Generalized Beam Functions 

In SCET the PDFs are defined as matrix elements of collinear fields with collinear proton 
states and are equivalent to the PDFs in QCD [43]. We illustrate this by giving the definition 
of the quark PDF in SCET and QCD respectively, 

jf CET (fc-/ P -,/,) = 0(fc-)(p»(p-)|x„(o)<y(fc- -Vn) |xr»(o)|p»(p-)) , (3.1) 
^(k- lp~^) = e[k-) ^e^-y + i\ Vn {p-)\^ 

The Fourier transform over the light-cone separation y + of the two fields in QCD is replaced 
by a delta function of the conjugate large label momentum in SCET. The Wilson line in the 
QCD definition is absorbed into the SCET fields, see eq. (2.1). In SCET the PDFs do not 
contain full QCD fields but collinear fields, which all belong to the same collinear sector. Each 
collinear sector is a boosted copy of QCD. Collinear fields may have zero-bin subtractions, 
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but these are not present for the PDF, since it is insensitive to soft radiation. (QCD fields 
have no such subtraction.) 

For {t, k 2 } 3> Aq CD , BFs can be related to the PDFs by performing an operator product 
expansion in Aq CD /t and Aq CD /fc^, 



BJs,x,k]_,n B ) 



B^(t,X,k ± ,HB) 



j 

£ 

i 
+ 



dx' 



x 

1 dx' 
x' 



T z <ii (*> 4' £i> ^ B ) h ( x '> Vb) 



, , )( A QCD A QCD \ 
" ' i k 2 ' 



(3.2) 



d- 2 



k 2 + 2 



9j \ x' ' ±, ^ B 



fj(x',fJ, B ) 



1 + 



A? 



A* 



QCD iv QCD\ 



t 



k\ 



where the index j sums over parton flavors. As stated in the introduction, we have assumed 
that t and k 2 are parametrically of the same size and thus get "integrated out" at the same 
time. The matching coefficients Xij and X%j will contain logarithms of k±/t, which are small 
because of this assumption. 

We determine the matching coefficients in eq. (3.2) by replacing the proton states in the 
BFs and PDFs by a collinear quark or a collinear gluon. At one loop, where at most one 
emission takes place, the invariant mass of the ISR is r 2 = 0. Eq. (1.1) then becomes an 
equality, 

k 2 = {1 ~ x)t . (3.3) 

x 

It is now straightforward to obtain the one-loop matching coefficients lij from the calculation 
of the one-loop standard beam function in refs. [14, 25], by including the extra perpendicular- 
momentum 5-function in eq. (2.5) and using eq. (3.3). The one-loop calculation of the gluon 
BF with uncontracted Lorentz indices is given in app. B, from which we obtain X^. The 
tree-level and one-loop matching coefficients are given by 



T§\t,x,kl,n B ) = -5 ij 5(t)8(kt)6(l-x) 



1 



7T 



i§\t,x,kln B ) = o. 



93 

X\l\t,x,kl,ti B ) 



I^\t,x,kl,j2 B ) 



a s {pL B )Ci,, 



2ir 2 



G(x) 



W 1 



5(l-x)5(kl) + \l 



k Vb 



Pu(x) 



x5 k 



2 (l~x)t 



+ S(t) 5{k\ 



TT 



a s {^B)Cij 
2tt 2 



6(x) 



2 -£o 



5 k 



Pu{x)\nx - '—5{l - x) +Zj l {x) 
o 

(l-x)t 

X 



+ 5(t)5(kl)ln 



1 — x 



P^ + S^dikDZtAx)), 



i^j(t,x,kl,i2 B ) 



a s (n B )C, 



y.i 



7T^ 



6{x)6{l-x 



1 — X 



(3.4) 
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where the plus distributions £ n are defined as 
C n (x) 



~6(x) ln n x~ 


= lim 


~e{x 


- (3) ln n x 


X 






X 



+ 8(x - (3) 



ln n+1 /3l 
re + 1 



(3.5) 



and satisfy the boundary condition L dx C n (x) = 0. The expressions in eq. (3.4) involve the 
color factors 



Cqq — Cgq — Cp 



a 



99 



19 



(3.6) 



the splitting functions 
P qq (x) = (1 + x 2 ) £ (1 - x) , 
P gg (x) = 2xC (l-x) + 2 9(l-x) 



rl 



P qg (x) = 6(1 - x) [x 2 + (1 - x) 2 ] 
l + (l-x) 2 



+ X(l-X) , Pgg(x) = 9(1 - X) 



(3.7) 



[where P qq is unconventionally defined without the 3/2 5(1— x) term to simplify eq. (3.4)] and 



Z 5 qq (x) = (1 + x 2 )d(l -x) + 9(1 -x)(l-x), 



^■qq( x ) 



2(1 - x + x 



2\2 



I 5 qg (x) = 29(1 - x) x(l - x) , 
l s gq (x) = 9(l-x)x. 



(3.8) 



Note that the cancellation of IR divergences in the one-loop matching of the standard BFs [25], 
immediately carries over to the generalized BFs as well. 

3.2 Comparison with Matching Coefficients of iBFs 

The conversion of our results to impact-parameter space is needed to compare with the 
matching coefficients for iBFs in refs. [10, 11]. Here we give the necessary equations both 
for CDR and CDR2 schemes. We will show that for our calculation CDR and CDR2 yield 
the same result in impact-parameter space (even though CDR is not suitable for momentum 
space, see sec. 2.2.2). 

For the two Lorentz structures in eq. (2.10) the transformation to the impact parameter 
y± in CDR is given by 



/ 



d-2 



9? 



k 2 d-2 



1 



( y±y± , jL 

\ y 2 ^ d-2)(l-e)(2-e 



d-2 
5 d - 2 (k ± -V n± ) 

rt2<f)2 



qFi 1 



X 



y±K 



2t52 



0F1 3 - e, 



(3.9) 



with d = 4 — 2e. The corresponding expressions for CDR2 follow directly by taking all 
perpendicular quantities in two dimensions and setting e = in these equations [except for 
the <7j_" /(d — 2) on the first line, which is just an overall factor]. The conversion needed for 
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the quark beam function is the same as that of the structure. To show how these results 
were obtained, we present the derivation of the first line of eq. (3.9): 



d d ~ 2 k ± 5 d - 2 (k ± - V n ±) 



I' 



d ~ 2 k\ e^ 1 



r(f) 



Ifci 



|4— rf 



[d - 2)ir d / 2 - 1 



(d - 2)7^/2-! ' 
d d - 3 ndkl e i|#xl|fex|cos0^2 _ £2j 



FA I 



0^1 



fin 



The angular integral was performed using 

2vr (d-3)/2 



d d ~ 3 n 



sin 



d-4 



(3.10) 



(3.11) 



where <j) is the angle between k± and y±. 

In sec. 2.2 we showed that the renormalization of the BF does not affect k±, so the UV 
divergences are multiplied by 8(k±). Since the PDF does not depend on k± (by definition), 
the IR divergences in the BF are multiplied by 8(k±) as well. Therefore any 0(e) pieces that 
may appear in eq. (3.9) due to the difference between CDR and CDR2 never get enhanced 
by 1/e divergences, 



Fi(l-e,0) = 1 =0^1(1,0) 



FA I 



FA I 



0^1 



+ 0(e). (3.12) 



Using eq. (3.9), it is now straightforward to compare our results with the matching 
coefficients for the iBFs in refs. [10, 11]. We agree with I gg , X qq and the matching coefficient 
for the new k± -dependent Lorentz structure in eq. (3.4). We disagree with I gq and I qg and 
discuss what we believe is the problem in refs. [10, 11]. Their calculation relies on the fact 
that in pure CDR the partonic fi/j is 



1 a s (y)Cij 
e 2vr 



e(x)p ij ( x ), 



implying that lij is simply the finite part of the one-loop beam function, 



B 



(i) 

i/j 



finite 



z, 



(o) 



f (1) + T (1) 



f 



(o) 



3/3 



finite 



(3.13) 



(3.14) 



However, the definition of the PDF, which yields eq. (3.13), requires one to average over the 
spins of the incoming parton. They do not do this correctly for B g j q and B q / g , where they 
sum over the spins and then divide by d — 2 (rather than 2) for B g j q and 2 (rather than d — 2) 
for B q /g. This incorrect averaging leads to a modification of the matching in eq. (3.14). For 
example, for B q / g 



B 



(i) 
q/g 



finite 



'q/9 



(0)' 

'9/9 



finite 



T (i) , a s (^)T F 



9(x)P qg (x). (3.15) 
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Including this additional term, their result agrees with our calculation. Note that since in 
ref. [25] and in app. A of the present paper the PDFs are explicitly calculated, these derivations 
are not prone to this problem. A similar issue occurs for I gq , but there in addition the sign 
of the 6(t) ln[(l — x)/x]P gq (x) term in X gq does not agree with the calculations in ref. [14] 
and our app. B. We were unable to identify the origin of what we believe is a sign error in 
Ref. [11]. 

3.3 Fully-Unintegrated Fragmentation Functions 

The FFs in SCET, which are defined in terms of collinear fields and collinear intermediate 
states \Xh) [26, 27], are equivalent to the standard QCD FFs that appear in factorization 
theorems at leading power, since the sum over the intermediate states is dominated by jet- 
like configurations [44]. As for PDFs, the zero-bin contributions cancel against each other in 
the sum over the graphs at each order, because these functions are insensitive to the scale 
associated with the ultra-soft radiation accompanying the parton [27]. (This is not the case 
for FJFs, where the zero-bin contributions do not completely cancel each other.) 

The matching of the FJFs onto standard FFs corresponds to an operator product expan- 
sion analogous to the case of the BFs. Our calculation at one loop in ref. [27] represents an 
explicit check that the additional measurement on the jet invariant mass does not spoil the 
cancellation of the IR divergences between FJFs and FFs. 

Starting from 

• Jz z ^ Z / 

the matching is performed by replacing the observed hadron in the FJF and in the FFs by 
a collinear quark or a collinear gluon. With at most one real emission, the perpendicular 
momentum pj_ of this parton is completely fixed in terms of s and z 

fl = 2(1 - z)s (3.17) 

which follows from momentum conservation and the on-shell condition p 2 = 0. As for the 
BFs, we can now straightforwardly obtain the matching coefficients from the calculation of 
the one-loop standard fragmenting jet function in ref. [27], by adding the delta function for the 
perpendicular momentum. There the quark fragmenting jet function was calculated using a 
gluon mass m as IR regulator. The gluon mass would in principle modify eq. (3.17). However, 
at the end of the calculation the limit m — > is taken to isolate the IR divergences. It is 
safe to take this limit ahead of time in eq. (3.17) since it does not modify the result for these 
divergences. 



1 + (^QCD Aqcd^ 



(3.16) 
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We find that the one-loop matching coefficients are given by 




x 6\pg ± - 2(1 - + <J( 5 ) P«(z) lnz - -5(1 - z) + 



2(2^) 3 




+ «5( 8 ) 6(j>g ± ) In (2(1 - 2)) P 3l (z) + «5( S ) 5(^)4(2) 



(3.18) 



where Cij, Pij(z) and lfj(z) were given in eqs. (3.6), (3.7) and (3.8). 
4 Examples of Factorization Theorems 

This section is devoted to illustrative examples of leading-power factorization theorems in- 
volving generalized BFs and FJFs. The first one concerns the transverse momentum Q±- 
distribution of a Drell-Yan lepton pair at perturbative values of Q± with a veto on hard cen- 
tral jets. This factorization formula can be straightforwardly extended to the case of Higgs 
production through gluon fusion (gg —> H) and Higgs Strahlung qq — > HV . The factoriza- 
tion theorem for gg — >■ H was used in sec. 2.2.4 to discuss the renormalization of BFs. The 
purpose of this paper is to describe the general picture associated with BFs, which suggests 
that they will play a role in describing transverse momentum distributions for events with 
any number of energetic well-separated jets, e.g. pp — > W/Z + n jets, once open factorization 
issues concerning for example the contribution from Glauber modes are settled. 

As a paradigm of factorization theorems involving FJFs, we will discuss the p± -distribution 
of the observed hadron in e + e~ — > Xh, where a cut on the thrust event shape is applied to 
restrict to dijet final-state configurations. Such a cut on thrust is utilized by the Belle collab- 
oration to study light-quark fragmentation in their on-resonance data, in order to remove the 
large 6-quark background. Here we consider the case of spin-averaged fragmentation. When 
polarization and azimuthal correlations between hadrons in back-to-back jets are taken into 
account, our framework can be applied to study Belle data on polarized light di-hadron frag- 
mentation in opposite hemispheres e + e~ — > dijets + 2h. 

4.1 Transverse Momentum Distribution in pp — > £ + £~ + jets 

We start our discussion by considering Drell-Yan (DY) production, pp — > £ + £~~, of lepton 
pairs with a large invariant mass Q. Factorization at leading power in Aqcd/Q has been 
discussed by Collins, Soper and Sterman, for any value of the transverse momentum of the 
lepton pair Q± with respect to the beam axis, namely both cases Q± ~ Q and Q± <C Q [45]. 
Since Q± equals the sum of the transverse momenta of the colliding partons entering the hard 
subprocess, the DY cross section is directly sensitive to the partonic transverse momenta. 



- 19 - 



The rigorous proof of the cancellation of contributions from Glaubers modes, which 
amount to initial-state hadron-hadron interactions taking place before the annihilation, that 
would spoil factorization [46], has been achieved in ref. [47] for the case of inclusive DY. 
Ref. [7] discusses an extension to the isolated DY process, pp — > £ + £~ +0 jets, where a central 
jet veto is imposed, restricting energetic ISR to be close to the beam axis. In this framework, 
the colliding partons are far from threshold (as in the inclusive case) and the collinear radi- 
ation is described by beam functions. Here we consider a more differential case, combining 
this with the study of transverse momentum distributions. We impose the central jet veto 
through a cut on the beam thrust event shape 7c m > defined in the hadronic center-of-mass 
frame as [14] 

r cm = J2\Pk\^ Vk] - (4.i) 

k 

The sum on k runs over all particles in the final state, except for the signal leptons, where 
pt and rjk denote the transverse momentum and rapidity of the particle with respect to the 
beam axis. The jet veto, T cra < 7^f -C Q, leads to large logarithms a™ ln m T^/Q (with 
m < 2n) in the cross section. 

If Aqcd Q± — V%mQ <C Q, we can write the following leading-power factorization 
formula following ref. [7] 

dc sr^ tt 1^2 \ f ixr f i. i. nii ( I- e~ Y t a +e Y tb 



dTcm dQ[ 



aoJ^HijiQ 2 ^) JdYjdt a dt b S%(T cm - e ~ ta + 6 |g Jdk a \dk b W 

6 [Ql-{k? ± +2\k±\\k£\ cos 0+fc 6 i)] Bi(t a , z a , ka±>») B 3^ z b , k b ±,v) , (4.2) 



where do is the Born cross section and <fi is the angle between and kir. The momentum 
fractions are given by 

Q Y Q -Y 



e , x b = — e , (4.3) 

"cm "cm 

with E cm being the center-of-mass energy and Y the total rapidity of the leptons. The 
sum in eq. (4.2) extends over the various quark flavors without mixed terms, i.e. ij = 
{uu, uu, dd, . . .}. 

The hard function H q q encodes the virtual effects at the hard scale Q due to the un- 
derlying hard process. It only depends on the large momentum components and not on the 
perpendicular momenta. The ultra-soft radiation is described by the soft function S q ^ , and 
the contribution of the soft radiation to Q± is power suppressed, due to our hierarchy of 
scales. Therefore H q q and the soft function Sg are the same as in ref. [7]. Only the BFs 
account for the recoil of the energetic initial-state radiation against the final-state leptons. 

Eq. (4.2) and the one-loop matching coefficients in sec. 3 provide all the ingredients 
needed to analyze up to NNLL the DY cross section in terms of 7cm and Q±, which can be 
tested for example against LHC data where the effect of resummation of large logs is expected 
to be important due to the large separation of scales. Because of our assumption ~ T C m Q-, 
the logarithms of T^ m Q/Qf in the cross section are not large. These are generated by the 
corresponding logarithms of t/k"j_ in the beam functions, through the convolutions in eq. (4.2). 
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Our approach provides two handles on the ISR, by simultaneously measuring its trans- 
verse momentum and beam thrust. This can be useful in the discovery and interpretation of 
new physics according to the authors of ref. [48]. For the example of SUSY di-squark (and 
di-gluino) production at threshold, their method leads to the determination of the sparticle 
mass. This entails a boost along the (beam) z-axis to the frame where the net p z of the 
final state radiation is zero. They determine the boost parameter from measurements on 
the final-state radiation, but this may also be achieved using beam thrust [8], which has the 
advantage of being unaffected by the presence of invisible final-state particles. In addition, 
their approach relies on the determination of the transverse momentum of the ISR. Our paper 
provides the tools for a reliable calculation of the cross section differential in beam thrust and 
the transverse momentum of the ISR, which is the input of their method. 

4.2 Measuring the Full Hadron Momentum in e + e _ — > dijet + h 

In ref. [27], the spin-averaged fragmentation of a light hadron was studied in the process 
e + e~ —> dijet + h. Experimentally [29], the dijet limit is imposed by a cut on the thrust event 
shape [49] 

T = max £ ^L^T> (4-4) 
l^i \Pi\ 

which we will use as well. In eq. (4.4), the sum runs over all particles in the final state and 
the value of t that maximizes T is called the thrust axis, which in the two-jet limit is along 
the direction of the jets. We will use the variable r = 1 — T, where r close to corresponds to 
configurations with two narrow, pencil-like, back-to-back jets. The maximum value of r = 1/2 
instead corresponds to a spherically symmetric event. The cut on thrust, r < r cut <C 1, leads 
to double logarithms a™ ln m r (m < 2n) in the cross section, that were resummed up to NNLL 
in ref. [27]. Here, we consider the situation where, in addition, the transverse momentum of 
the hadron pj^ with respect to the thrust axis is measured. For intermediate values of Ph±, 
namely for Aqcd "C Ph± ~ VtE cu1 <C E cm (in accordance with the SCETi power counting) 
the following factorization formula holds 

d(J v-^ Oq ^ < - ~ (~ s a + 



drdzdpfi ± ^ (47r) 2 



H(Q 2 ,fj,) Jds a ds b Q S r (Qt - ScL Q 2 Sb i M 



X [Oq(s a ,Z,Ph ± ,n) Jg(s b ,n) + J q (s a ,fl)g^(s b ,Z,p^j_,n) , (4.5) 



where <tq is the Born cross section. Since the hadron is produced by the energetic radiation 
inside a jet, assuming cancellation of Glauber contributions, the only modification required to 
the right-hand side of the factorization theorem in ref. [27] is the replacement of the standard 
FJFs by the generalized FJFs. The hard function H encodes virtual effects from the hard 
process at the hard scale E cm , the soft function S T describes the contribution to thrust due 
to ultra-soft emissions, and the jet function J is associated with the energetic radiation in 
the jet which does not contain the detected hadron. All these functions are the same as in 
ref. [27]. Together with the matching coefficients in sec. 3 we provide all the ingredients to 
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determine the cross section in eq. (4.5) at NNLL order. Since p^± ~ the logarithms of 

V^Q/Phi i n the cross section are not large. 

5 Conclusions 

In this paper we defined and studied the generalized BFs and FJFs. The former describe the 
distribution in the full four momentum k^ of a colliding parton taken out of the incoming 
hadron, for perturbative values of k^_. We considered t = —k + k~ ~ k^_ which avoids intro- 
ducing an additional hierarchy between these two scales. Generalized BFs are relevant for 
factorization theorems for exclusive iV-jet processes, when the measurement is sensitive to 
the transverse momentum of the ISR. We show explicitly that a proper definition of general- 
ized BFs requires ultra-soft zero-bin subtractions. We have argued that the renormalization 
of the generalized BF is the same as that of the jet function, and discussed an important 
subtlety in defining and calculating perpendicular momentum dependence using dimensional 
regularization. We have calculated the one-loop matching of BFs onto standard PDFs and 
corrected the non-diagonal coefficients X gq and X qg previously worked out in ref. [11]. 

We similarly discussed the generalized FJFs, which describe the full momentum de- 
pendence of a hadron h fragmenting from a parton with virtuality s. The perturbative regime 
s ~ pfij_ S> Aq CD is assumed. An example for which we give the factorization theorem in- 
volving generalized FJFs is e + e~ — > dijet + h, where a cut on the thrust event shape is used 
to impose the dijet limit and the full hadron momentum is measured. The Belle collabora- 
tion employs such a cut on thrust to remove £?-meson events from the data sample on the 
T(4S") resonance and study light-quark fragmentation [29]. The extension of our framework 
to include spin correlations, would allow one, for example, to study the Collins effect [50] . We 
showed that the renormalization of the FJF is equal to that of the jet function and obtained 
the one-loop Wilson coefficients for matching FJFs onto standard fragmentation functions. 
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Figure 3. Feynman graphs contributing to the quark BF at one loop. Graphs (b) and (c) have 
symmetric counterparts that are not shown but are included in their computation. Graph (e) vanishes 
in the Feynman gauge. 



A One-Loop Quark BF with IR and Rapidity Regulators 

We present detailed results for the quark generalized BF calculation at one loop, showing 
the structure of the IR and UV divergences, the zero-bin contribution and the cancellation 
of the rapidity divergences. To ensure an unambiguous calculation, we regulate the IR, 
UV and small rapidity region in the integral with distinct regulators. We use dimensional 
regularization for UV divergences, for IR divergences we use a gluon mass (or a quark mass 
for the mixing contribution) and to regulate rapidity divergences we use the 5-regulator [51]. 
Rapidity divergences arise in a loop- integral with eikonal propagators when k + k~ is fixed but 
either k + or k~ goes to zero, see for example [52]. The (5-regulator of ref. [51] takes care of 
such divergences, however, historically it has been thought of as an IR regulator. In ref. [27] 
we calculated the fragmenting jet functions and fragmentation functions employing this same 
choice of regulators. We find that the UV divergences extracted from this calculation of the 
generalized BF are the same as those of the jet function and the standard beam function [7, 25] . 
We will show that only after an ultra-soft zero-bin subtraction the rapidity divergences are 
cancelled in the BF and the IR divergences match with those of the PDF, as expected. 
Consequently we find that the (quark) BF, without zero-bin subtractions, contains rapidity 
divergences. 

The graphs contributing to the quark BF at one loop are shown in fig. 3. Including the 
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symmetric counterparts for (b) and (c), the result for each graph is: 
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Graphs (b) and (c) have also non- vanishing zero-bin contributions: 
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Graph (e) vanishes in Feynman gauge. The sum of these diagrams yields B q / q , 
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where the first line of the second equality contains the UV divergences, the second line the IR 
divergences that match with those of the PDFs (which are given below) and the last two lines 
show the finite contribution from BFs to the matching. As expected, the In 5 terms cancelled 
out after including zero-bin subtractions. 



- 24 - 



Graph (f) contributes to the mixing term in the matching 5 , 
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with the splitting function P qg given in eq. (3.7). The crossed version of graph (f) does not 
contribute for positive values of x. 

If the ultra-soft zero-bin subtractions are not performed, the ln<5 will remain in the 
final answer, indicating the presence of rapidity divergences. If they are interpreted as IR 
divergences (as is typically done) then the BFs cannot be matched onto PDFs, since the IR 
divergences of the BFs do not match up with those of the PDF. Alternatively, if these are 
not considered IR divergences and absorbed into the BF's renormalization factor (as in a 
pure dimensional regularization calculation) the matching onto PDFs is still feasible. In this 
approach there is a lack of distinction between rapidity and UV divergences. However, this 
alternative does not apply to our BFs, since there is a definite measurement that restricts the 
rapidity and prohibits rapidity divergences. 

In order to perform the matching, we need also the PDFs calculated with m and 5 
regulators. Our final result is 
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(A.6) 



where the (^-dependence cancels. Here m denotes the gluon mass in f q / q and the quark mass 
in fq/g as was the case for beam function. We obtain the same matching coefficients as in 
sec. 3. 



B One-Loop Gluon BF in Momentum Space 

In this section we calculate the gluon BF with uncontracted indices fj,, u, using CDR2. The 
diagrams are shown in fig. 4. We work in pure dimensional regularization, which sets the 
virtual diagram (c) and the wave function diagram (d) to zero because they contain no 
dimensionful quantity. We will use Feynman gauge, which makes diagram (e) vanish because 
the gluons emitted by the Wilson lines have n-polarization and n ■ h = 0. By only summing 
over the physical polarizations of the final state gluon [see eq. (B.3)], diagram (b) vanishes as 
well. Thus only diagram (a) contributes to B 9 g and only (f) contributes to B g . 

5 Here m is quark mass rather than the gluon mass, but this neither affects the renormalization nor the 
matching. 
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Figure 4. Feynman diagrams contributing to gluon BF at one-loop. Diagrams (b) and (c) have 
symmetric counterparts that are not shown. Diagram (d) represents self-energy corrections. Diagram 
(b) vanishes because we restrict the polarization sum to physical polarizations. Diagrams (c) and (d) 
vanish in pure dimensional regularization and diagram (e) vanishes in the Feynman gauge. 



We will start with diagram fig. 4(a), first giving the expressions and then we will comment 
on the calculation and several of its subtleties. 
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The plus distributions C n are denned in eq. (3.5) and the splitting function P gg in eq. (3.7). 
Following the standard gluon BF calculation in ref. [14], we will average over the colors and 
polarizations of the incoming gluons, 



pol 

For the intermediate gluons we restrict the polarization sum to physical polarizations, 

^ e A (£) £p (£) - -g Xp + — + — - . (B.3) 

pol V ' 

The polarization sum for the incoming gluons is simpler because = 0. 

In the first step we performed all the Lorentz contractions and integrated most of the 
delta functions. It is convenient to decompose £j_ = 1% +£e, where £2 lives in 2 dimensions and 
£ 6 in — 2e dimensions. In CDR2, the transverse momentum measurement does not constrain 
£ e and its integral produces the required divergences, 

W U{-£l){-£ 2 t y 1 -*. (B.4) 



(2tt)-* r(-e) 

The £e£ v t Lorentz structure produces a — £ 2 ge U /(2e), where g^ v is the — 2e-dimensional part 
of the metric tensor. 

To clearly separate the variables in the third expression, we swapped k 2 for u = xk?/[(l— 
x)t], which is kinematically restricted to < u < 1. We included the Jacobian for changing 
the variable of the distribution from k 2 to u, since distributions are tied to the corresponding 
measures, 

B(kl,...)dP ± = B(u,...)du. (B.5) 
We then use the distribution identity 

= - 1 5 (z) + C (z) - ed(z) + 0(e 2 ) , (B.6) 

to obtain distributions in t/fi 2 , l — x and 1—u. (The du/dk]_ should not be expanded.) Based 
on the one- loop kinematics we would expect the expression to be proportional to 5(1 — u), 
but the expansion also produces a Cq(1 — u) and C\(l — u) as well. However, all these terms 
are multiplied by a 5(t) or a 5(1 — x), implying that k 2 = 0. The situation is similar to having 



-27- 



a <fi dependence at = in spherical coordinates, where all choices of 4> are equivalent. We 
are therefore free to replace this u dependence by a delta function times its average value 
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After this, we change variables back from u to k]_, where the Jacobian du/dk]_ gets absorbed 
into (5(1 — u) to yield 5[k"j_ — (1 — x)t/x] in the final expression. 

We write the final result in terms of independent Lorentz structures. For the new tensor 
term we have 
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0, 
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because if k"j_ = there is no preferred direction and kj_k^_/k]_ — > —g^/2. This observation 
eliminates the divergent contribution for this Lorentz structure, which cannot be present 
because this tensor first appears at one loop. In the final result the expressions multiplying 
<7j_ are the same as for the standard BF [14, 25]. The structure, that appeared in 
intermediate expressions, is absorbed into this piece. This is not surprising because there is 
no preferred direction in the — 2e-dimensional part of space. For all intents and purposes, 
we could therefore have taken the fi and v indices in 2 dimensions from the beginning [the 
factor of 1/(1 — e) in /(l — e) will never affect the renormalization nor the matching]. We 
reproduce the known renormalization and matching for the g^ piece and obtain the matching 
coefficient for the £;j_-dependent Lorentz structure. 

The calculation of the diagram in fig. 4(f) with external quarks instead of gluons is very 
similar 
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This result also agrees with the matching obtained from the standard beam function calcu- 
lation. All the matching coefficients are collected in sec. 3.1. 
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